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Abstract. We consider some parametrized classes of multiple sums first studied by 
Euler. Identities between meromorphic functions of one or more variables in many cases 
account for reduction formulae for these sums. 



1 Introduction 

A somewhat unlikely-looking identity is 

oo n— 1 oo ^ 

y - y 1 = y - , (i) 

^— ' n(n — x) ^— ' m — x ^— 4 n 2 (n — x) ' 

ra=l v ' m=l n=l v ' 

valid for all complex x not a positive integer. For x — 0. (|Tj) becomes ((2, 1) = £(3) — the 
most central Euler sum identity — as discussed below. 

In this note we begin with an ab initio proof of and then explore various conse- 
quences and extensions. We conclude by studying other generating function identities of 
which 

n sr~y 2m 
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y = (coth(7rj/) + coth(27rj/)) V — =^ ^ (2) 

n 2 + 4y 2 v v v ^ (n 2 + y 2 )(ra 2 + 9y 2 ) KJ 

is a pretty example. 

Though we do not belabour the point, all our work was assisted by the use of computer 
algebra systems — and some of it would have been impossible, at least for us, without such 
tools. This is true of both our discoveries and of our proofs. The joys of such symbolic 
and numeric computation . . . and much more are discussed in detail in pQ and [2] . 



2 A Parametric Euler sum 

We begin by defining various special functions which we shall exploit in the sequel. For 
Re(x) > 0, the gamma function and its logarithmic derivative are defined by 
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T(x) := / e-H^dt, V(x) := 
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r(x) 
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Note that 

x ( N 1 \ 

E — r = — ^(T — x) — 7, where 7:= lim ( > logn I 

n(n-x) K ' N^oo\^n ) 

n=l v ' v n=l 7 

is Euler's constant. 

Recall that the classical Riemann zeta function is defined by 

n=l 

Correspondingly, 

1 x 1 n ^ 1 

n>m>0 n=l m=l 

defines a double Euler sum. The two-place function @ was first introduced by Euler, 
who noted the reflection formula 

C(s, t) + C(t, s) = C(s)C(t) - C(s + t), Re(s) > 1, Re(t) > 1, (4) 

and the reduction formula © below. An obvious extension of (J3J) is 

C(si, s 2 , . . . , s N ) := ^2 



Si So S/tf ' 



ni>n2>--->njv>0 

which defines an Euler sum of depth N and weight s«. Some authors reverse the order 
of the variables. 

Euler sums may be studied through a profusion of methods: combinatorial, analytic 
and algebraic. The reader is referred to |2J Ch. 3] for a concise overview of Euler sums 
and their applications. 

We now prove identity (JTJ) directly. 
Theorem 1 If x is any complex number not equal to a positive integer, then 

n-l 



y y 1 =y t 

^— ' n(n — x) ^— ' m — x ^— ' n 2 (n — x) 
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Proof. Fix x G C \ Z + . Let 5 denote the left hand side. By partial fractions, 
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Now for fixed m 6 Z + , 

N 



^— ' \n — m n J iv-^oo ^— ' \n — m n J ^— f n jv-^oo ^— ' jV — n + 1 

n=m+l x 7 n=m+l x 7 n=l n=l 

n=l 

since m is fixed. Therefore, we have 

oo m 1 oo n— 1 oo ✓ 71 1 n— 1 1 
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For x = v— T, equation Q becomes the pair of tangent sum evaluations 

oo n— 1 oo _. 

El v ^ mn — 1 x ^ 1 , — 
n(n 2 + 1) ^ m 2 + 1 = ^ n(n 2 + 1) = 7 + 6 ^ + 



and 
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El \ - n + m \ - 1 it t / 1 / — 7\ 
> = > = Im^fl + V— 1) 
n(n 2 + l)^m 2 + l ^-f n 2 (n 2 + 1) 6 V ; 
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7T 2 7T COtll 7T — 1 

~ y 2 ■ 

Setting x = in equation (JTJ) gives £(2, 1) = C(3). More generally, differentiating ((TJ) 
A; times with respect to x produces a corresponding formula for ((k + 3). For example, 
differentiating once and setting x = (equivalently, comparing coefficients of x on both 
sides) produces 

C(4) = C(3,1) + C(2,2). 
Since, by the reflection formula (jlj), £(2, 2) = [C 2 (2) — £(4)]/2, we also evaluate 

1 7T 4 

C(3,1) = 7 C(4) 



' 360 



This is the first case of a remarkable identity discovered by Zagier, and first proved in j3] . 
See (dU) below. 

Integrating equation (jlj produces the following corollary. 
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Corollary 2 For all complex x not equal to a positive integer, 

oo oo n-1 , , x 

z — ' n z V n/ ' n ^— ' n — m \l — x/nj 

n=l n=l m=l x ' 7 

For x = a/— 1, the imaginary part of (j3J) leads to 

I 1 arctan ( m ) = arctan(l/n) = (-1)" + 
^ n ^ m \n 2 — mn + lj ^ n 2 ^ 2n + 1 

n=l m=l x 7 n=l n=0 

where the last evaluation comes from writing arctan(l/n) = L n(x 2 + n 2 )~ l dx and ex- 
changing the order of summation and integration. 



3 Euler's Reduction Formula 

Euler's reduction formula is 

s-2 

C(s, 1) = -s((s + 1) - -J2 C{k + 1)C(* - 1 < s G Z, (6) 



fc=i 



which reduces the double Euler sum £(s, 1) to a sum of products of classical Riemann 
(^-values. Another marvellous fact is the sum formula 

C (ai + 2, a 2 + 1, • • • , a r + 1) = C(r + s + 1), (7) 

a;>0 

valid for all integers s > 0, r > 1. It is easy to show that (JJJ) is equivalent to 

r oo 

E -rUr^ = JI^—m ^ Z+ - ( 8 ) 

k 1 >k 2 >--->k r >0 j=l J n=l v ' 

The first three non-trivial cases of © are C(3) = C(2, 1), C(4) = C(3, 1) + C(2,2) and 
C(2,1,1) = C(4). 

The ordinary generating function of the sequence {((s) : 1 < s G Z} is 

oo oo 

Z{x) := COO x s - 1 = = " - 7- 

z — ' z — ' nm — x) 

s=2 n=l v 1 

One can check that Euler's reduction © is equivalent to the ordinary generating function 
identity 

x n— 1 ^ x j f ^ 1 f C 

nfn — x) n(n — a;) 2 2 | n 2 (n — s) 2 \ nfn — a:) J | ^ ^ 

n=l v y fc=l n=l v ' Kn=l v 7 v n=l v J/ ) 



This relies on observing that the right hand side of © involves the square and the 
derivative of Z{x). In turn, (j^J) is equivalent to (JJJ). This equivalence may be demonstrated 
as follows: Let 
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- <7„.= Y^ - .<7„.= V^ - 2 . 
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z — ' n(n — x) z — ' k[k — x) 
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and 

n-l 
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Combining these yields the desired equivalence. 
Correspondingly, equation ((TJ) is equivalent to 

C(s + 3) = C(2 + 1 + 6), < s G Z, 

a,6>0 

which is an inversion of Euler's reduction formula (jfijL and simultaneously recaptures the 
case r = 2 of equation (JJJ). 

Thus, Theorem [TJ has established all of these results directly, without using analytic 
methods, or sophisticated partial fraction decompositions, as is usual. The first case, once 
more, is £(2,1) = ((3). 



4 More Generating Functions 



Euler sums arise very naturally and it is perhaps a trick of time that C{2n) is viewed as 
more natural than 



n 

C({2} n ) := C(2^. ._^) = lb 



k 2 ' 

ki>k2>->k n >0 j=l i 



Indeed, Euler's infinite product for the sine function is precisely equivalent to 

= £(-!). f«2}. 
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n=0 



r 2n 
fnl i 



5 



from which we may deduce that 



n 2n 



Similarly, with uj a primitive cube root of unity, 



oo 00 / 3 \ 

!;<«»}.)*"• =n i+ f 

ra=0 n=l ^ 



r(i + x)r(i + ujx)v{\ + lu 2 x) 

1 + C(3)x 3 +Qc 2 (3)-^) * 6 + 0(* 9 ), 



and allows one to show that C({3}n) is always in the ring generated by the numbers ((3k) 
(fe = l,2,...,n). 

By various methods, one can show that 

C({3} n )=C({2,l}„) 

for all < n G Z, while a proof of 

C({2,l} n )=2 3 "C({2,l} n ), 0<neZ (10) 

remains elusive. The bar over the 2 on the right hand side of ffTTTf) signifies that in the 
summation, terms of the form 1/k 2 are to be multiplied by (— l) fc . Such sums are called 
alternating Euler sums. Only the first case of (fTUj) . namely 

00 ^ fc— 1 ^ oo ^ jhj, fc— 1 ^ 

E p? E - = 8 E -p- E - (= c( 3 )) 

fc=l m=l fc=l m=l 

has a self-contained proof E] ■ Indeed, the only other proven case 1 is 

00 ^ A;— 1 ^ m— 1 j p— 1 ^ oo / ^xj^ fc— 1 ^ m— 1 / ^ p p— 1 ^ 

EpE-E^E- = 64 E-^E-E— ^E- (=c(3,3)). 

fc=l m=l p=l r q=l * fc=l m=l p=l ^ <j=l ^ 

There has been abundant evidence amassed to support ()10|) since it was first conjec- 
tured |3] in 1996. For example, very recently Petr Lisonek checked the first n < 85 cases 
to 1000 decimal places in about 41 hours with only the expected roundoff error. And he 
checked n = 163 in ten hours. If true, (fTUj) would be the only known identification thus 
far of an infinite parametrized class of alternating Euler sums with a corresponding class 
of non-alternating Euler sums. An intriguing reformulation of (|10|) is stated below. 

1 This is an outcome of a complete set of equations for MZV's of depth four. 
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Conjecture 1 Define a sequence of polynomials a n = a n (t) for positive integers n by 
a\ = a 2 = t 3 and 

n(n + l) 2 a n+2 = n(2n + l)a n+1 + (n 3 + (-l) n+1 t 3 )a n , n>l. 

Then 



( t 3 \ 
lim a n = t 3 H(l + — ). 

n=l v ' 



A more recondite generating function identity is 

4n cosh(7ra;) — cos(7ra;) 



J]C({3,1},.) ' 

n=0 



7T 2 X 2 



which is equivalent to Zagier's conjecture (subsequently proved): 

C({3,1}»)= (^yy, < n e Z. (11) 

The proof of (fTTjl (see |2J p. 160], jl], 0) devolves from a remarkable factorization of the 
generating function in terms of Gaussian hypergeometric functions: 

oo 

C({3, l}n) x 4 " = 2 F l (t, -t; 1; 1) 2 F l (tt, -it; 1; 1), 

n=0 

where £ = (1 + i)x/2. 



5 Further extensions 

Relatedly and more centrally, Euler's reduction formula (JBJ) has many extensions. For 
instance, with even a > and odd 6 > 1 with a + b = 2N + 1, one has [3] 

C(a, b) = C(o)C(6) + ~{ f + ^ - 1 k(a + 6) 



2 1V a 

N-l 



'^{{a 2 -l) + { b 2 -l)}^ + 1 ^ + b - 1 - 2 ^ <12> 



r=l 

and hence by the reflection formula (j3J) 

ab,a) = - 1 -{l + ( a+ a b )}C(a + b) 

N-l 



(13) 

+ £{(„-l) + «* + !)«. + »- l-2r). 
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Although £(1) is undefined, Euler's reduction (0) implies that (|T2Jl also holds in the case 
b = 1 if we interpret £(1) = 0. Thus, (|12j) is indeed an extension of (JHJ). 

We next recast (|12|) using generating functions, and then develop some further con- 
sequences and interesting special cases. We initially keep b = 2t + 1 fixed and want to 
express 

oo oo oo 2 (iV-t-l) n ~ l 1 °° 1 n_1 1 

)= ^2 n 2(N-t) ^6 = X] n 2 _ x 2 ^6' 

JV=t+l JV=t+l n=l m=l n=l m=l 

We suppose that x 2 is not a positive integer and set D := — . The requisite component 

ax 

generating functions are 

oo oo oo 2(W-t-l) 00 J 

*■ ) = n 2(N-t) = n 2 _ x 2 ' 

iV=i+l 7V=i+l n=l n=l 



^2 C(a + b)x( } = n 2(N-t)+b = E n 6f n 2 _ 3.2 j 

iV=t+l JV=t+l n=l n=l v 

^ C(a + 6) ^ + & V2(iV-^l) =x -2 J2 C(2iV+l 
r— +_i_i \ / at— 



—x 2N+l 



9^,2 / / 



W 7 + 



2x 2 ^— ' \ (n — a;) b+1 (n + x) 

n=l x v ' v 

00 iV-1 / 9r \ 

J2 ^ N - l - 1] E U _ 1 c(2r + 1)c(2iV ~ 2r) 

/=t+l r=l ^ ' 



6+1 / x 2 



N=t+1 



1 \y Z^ 1 / X 2 

^ n 2 - x 2 ) ^ (b - 1)\ \n(n 2 - x 2 ) 

n=l ' n=l v / \ \ / 



00 ^ \ °° 1 

5j „2 _ ~2 ) ^2 9 



1/1 1 

+ 



' n 2 — x 2 J ^— ' 2 \ (n — a;) 6 fn + 

n=l ' n=l x ' v ' 



and 



£ ^-*-D E( a _ 1 C(2r+ 1)C(2JV - 2r) 

\r — +J_1 1 — :i \ / 

* 00 r,b-2m / 

m=l n=l v 7 v 

771=1 71=1 VV ' 



N=t+1 r=l 



n(n 2 — x 2 ) 
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6+l-2m ^ _|_ a; ^6+l-2m 
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m=l n=l N 



x) 2m {n + x) 2m i 

Combining these, with a fair amount of care, yields a generating function for fixed b: 



n — 1 

^— ' n z — x l ^— ' m° ^— ' 

n=l m=l s=l 



oo ^ ^ oo ^ 

C(2t + k 2 -x 2 ~ 2 ^ n 2m ( n 2_ x 2) 

=1 n=l v y 



oo 



1 n f 1 



i 

+ 



2 ^2x 2 l(n-x) 2 * +2 (n + x) 2t+2 n 2t + 2 



9 ( £ r,2 _ ~2 ) ' £ 



1 1 

+ 



2l^n 2 -xV ^ 1 (ra - x) 2 * +1 (n + x) 2t + l 

\ n= \ / n=l v v ' v ' 

t °° 1 f 1 1 

£ C( 2 * + 2 - 2m ) Yl 2x \ (n- x) 2m ~ (n + x) 2m 

m=l n=l v v 



Thus, 



1 71 COt (tTx) f 1 1 



^2 _ ™2 /_^t mb d-r. ^—^ 



+ 



n*-x'^m» Ax ^ \( n - x ) 2t+1 ( n + x ) 2t+1 n 2t + l 

n=l m=l n=l v v / \ / 

1 » n -2t-l 1 ^ - f ;| | 



using 

7rcot(7ra;) 1 1 1 



n=l m=l n=l v v ' v ' 

1 °° 1 1 

2x 2f"^^ and ^ = "2- 

fc=i 

Now we sum over the odd parameter, b, and obtain a two- variable ordinary generating 
function expressible as: 

n-l 



ny(4y 2 — x 2 ) / (n 2 — y 2 ) 



n=l m=l a s>0,t>0 



2 °° C ° t{nx) ^ '((n + a;) 2 - t/ 2 ) ((n - x) 2 - y 2 ) 



n 
oo 



+ — xcot(7ra;) 



((n + x) 2 — y 2 ) ((n — x) 2 — y 2 ) 
^ ((n + re) 2 - y 2 ) ((n - re) 2 - y 2 ) 



+ ^ycot{7ry)Y^ 
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-ly ^ . (14) 

2 ^ (n 2 - y 2 ) (n 2 - x 2 ) K ' 

Attractive specializations come with x = ±y, x = ±2iy, x = and, after dividing by y, for 
y = 0. For example, at (1/2, 1/2) where the right hand side has a removable discontinuity, 
we deduce that 

oo n— 1 o 

El \ ^ m 7r 1 

4n 2 - 1 ^ 4m 2 - 1 ~ 64 ~ 16' 

n=l m=l 

Note that the second, third and fourth terms on the right hand side of (I14|) are 
expressible as 

7T . .^(l+x + y)-^f(l + x-y)-^f(l-x + y) + ^f(l-x-y) 

L 2 := — — xcot(7ix) , 

4 Axy 

7T . . *(1 + x + y) - ^(1 + x - y) - ^(1 - x + y) —x - y) 

^3 ■= tzV cot{7Ty) 



£ 4 : 



2 4xy 
1 - {*(1 + x) - -y) + *(ar) - *(y) + vrcot(vrx)} y 



4(y 2 — x 2 ) 

respectively while the first term has a corresponding evaluation. It is 

2 . . *(l-x-y) + *(l + x-y) + #(l + x + y) + *(l-x + y) 



Si := y 7r cot (/T 

— y7T COt(7Tx) 

— y 2 n cot(7rx) 



8x(x 2 — Ay 2 ) 

^(1 + x + y) - ^(1 + x - y) - ^(1 - x + y) + ^(1 - x - y) 

16(x 2 -4y 2 ) 

*(2/ + 1) + *(1 - 2/) 



4x(x 2 — 4 y 2 ) 

We finish by expressing a form of this generating function concisely as: 
Theorem 3 For all x and y with squares not equal to negative integers, 
n ^-k 2m \ 1 



\ n 2 + y 2 r// 2 ■ 7-' / //- - ./- 2 

n=l N m=l ' 



x 



2 . > \ ^ ri 



— 4y 2 )irx coth(7rx) 



^ (n 2 + y 2 ) {(n 2 — x 2 + y 2 ) 2 + (2nx) 2 } 



+ (27ry coth(/Ty) + tcx coth(7rx)) 



' (n 2 — x 2 + y 2 ) 2 + (2nx) 2 

Proof. Replace x by ix and y by iy in fll4j) . and then factor denominators and regroup 
the terms as needed. QED 
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Letting x and y approach zero yields C(2, 1) = £(3) again. Setting x = 2y pro- 
duces (J2J) — the identity with which we began. 

To conclude, we challenge the reader to explicitly obtain the corresponding two variable 
generating function for ((2t + 1, 2s). 
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